233 


SH 2.1-9 

DIFFUSION IN THE CHROMOSPHERE , AND THE COMPOSITION OF THE SOLAR CORONA 

AND ENERGETIC PARTICLES 

Sylvie VAUCLAIR 

Observatoire du Pic-du-Midi et de Toulouse, France 
Jean-Paul MEYER 

Service d'Astrophysique, Centre d' Etudes Nucleaires de Saclay, France 


1. Introduction. Composition observations, (i) m the Solar Photosphere, 
and (li) m the upper Transition Region (TR) and Corona (as well as in 
Solar Wind, SW, and Solar Energetic Particles, SEP) imply a change of com- 
position of the solar atmosphere somewhere between the Photosphere and the 
upper TR (T >70000 K) /7, 8, 25, 26, 46/. Heavy elements with First Ionization 
Potential (FIP) > 9 eV (hereafter "high-FIP elements") are ~4 times less 
abundant in the TR and Corona than in the Photosphere, as compared to both 
Hydrogen and heavy elements with lower FIP ("low-FIP elements"). 

These observations suggest a separation between neutral and ionized 
elements in a region where the high-FIP elements are mostly neutral, and 
the low-FIP elements ionized /25,26/. This occurs in the Chromosphere at 
altitudes above 600 km and below 2000 km above Photosphere (in the average 
quiet chromospheric model VAL C of /43/jUsed throughout this paper; fig. in 
/ 1 3/ ) . Here we shall investigate whether diffusion processes can explain 
the observed change in composition. 

2. The diffusion equation - Possible processes. When a gas mixture is sub- 
ject to pressure, temperature or concentration gradients, or to selective 
forces, the gas components diffuse with respect to one another. This process, 
believed to be the reason for chemically peculiar A and B stars /e.g.28, 

29, 27, 45/, has already been invoked to explain the He depletion in the Solar 
Wind /19, 27, 12, 13/. Thermal diffusion in the TR has also been studied by 

/9, 30, 42, 36, 33, 34, 35, 12, 13/. The diffusion velocity for a gas with negligible 
abundance in a gas mixture may be written /6,10/ 

v = D ( Vln c - 1 l y In p - k T v In T + F/kT ) 
where D is the diffusion coefficient, k the pressure diffusion factor (the 
V In p term is equivalent to the gravitational settling) , k^ the thermal dif- 
fusion factor, F any selective forces (e.g. the radiative force), and c the 
concentration of the element. 

In solar-type stars, the radiative acceleration is never large enough 
to compete with gravitational settling /28/ : the last term is negligible. 

In the chromospheric plateau, the temperature gradient is small, and the 
thermal term is also negligible. In the TR this term is large. However, 

all elements become 100 % ionized as T increases and the neutral-ionized 
element separation therefore cannot work there. As the concentration grad- 
ient is also generally small the diffusion velocity reduces to 
v = D k p V In p (gravitational settling term) 

In case of partial ionization of an element, the atoms spend part of 
their time in the various ionization states, and the resulting velocity may 
be written, in first approximation /29,1/, as v D = »?k v Dk » where t\ k is the 
fraction of atoms in the ionization state k, and v D ^ the corresponding dif- 
fusion velocity. 

If the elements have to diffuse across a magnetic field, the tansverse 
diffusion velocities of the ions are reduced by a factor of (1+u t col ) _ , 
where w is the Larmor pulsation and t col the characteristic time between two 
collisions /44, 1/. If this factor is large, the transverse diffusion of ions 
is stopped, while the diffusion velocity of neutrals is unchanged. Thus the 
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diffusion velocity of any element reduces to v D = v n v Dn , where the~ 
index n refers to the neutral fraction. 

Now# a large fraction of the solar surface is always covered with 
horizontal magnetic fields anchored in the network, from altitudes of ~700 
km /3,11,47, 14,20,1s/ up to above the Ha emitting region (>1500 km)/22,23/. 
Their typical intensity is ~ 5 to 25 G. At a typical chromospheric tempera- 
ture of ~ 6000 K, the factor t CQ ^ « (4.6- lO^/njj) 2 for ions becomes >1 
for altitudes above ~ 950 km , reaching values of 10 and 100 at /v 1100 
and ~1300 km. Above these altitudes, ions are effectively prevented from 
settling gravitationally, while neutrals are not. 

3. The diffusi on coefficient for neutral heavies. The diffusion coefficient 
D for neutral heavy nuclei of mass A on the dominant neutral H gas of the 
chromosphere can be expressed as / 6 , 10 / 

D = 3\Z?ir / 1 6 * 1/ (njj • c e ff ) . -\/kT/m*" 

D = 4 . 287 • 10 • 1/ (njj *cr ^5 )*\/T»\/l + l /A (cm 2 s“l) 

where n H is the neutral H number density, T the temperature, m*= mgA/(A+l) 
~mfj the effective mass of the interaction, a eff an "effective" collision 
cross section to be discussed below, and < 7^5 = ^eff in units of 10 “^ cm 2 . 
(The diffusion coefficient in ionized H is not very different; a Small ad- 
mixture of H therefore does not affect the picture) . The effective cross- 
section <r eff can conveniently be expressed as <r Qff = tt r* 2 ^ 1 ' 1 '* , where 
r* is some typical range of the interaction potential, and q(1,1)* 1s 
reduced collision integral describing the collision kinematics in the actual 
potential /6 , 10,32,41 ,31/. 

The noble gas-H interaction potentials are characterized by a steep 
repulsive part and a tiny Van der Waals attractive well, whose depth e ~ 10 -2 
eV is extremely small as compared to the mean thermal energy kT = 0.52 eV of 
the atoms in the ~6000 K chromospheric gas (fig.l). So, only the repulsive 
part is important. The key information on these potentials -comes from /5/, 
whose results have been by and large confirmed by more recent work / 4 , 39 , 
40,17,38/. The repulsive part of the potential is most realistically fitted 
by a modified Buckingham- type potential /21,5,6,10/, and we have adopted the 
potential parameters of /5/. On this basis the collision integrals fid/ 1 )* 
have been interpolated from the tables of / 21 / (reproduced in / 6 , 10 /), yiel- 
ding the effective cross-sections <r eff given in Table 1 (or the effective 
collision radii r e ££ such that ^* e ff = w r eff^ shown in fig.l). 

In reactive element-H potentials, the tiny Van der Waals minimum is 
completely swamped out by a~4 eV attractive well associated with the cova- 
lent bond (fig.l). The characteristics of these potentials have been derived 
from the data tabulated in /18/, using the formalism of /16/. They can be 
reasonably well fitted with Morse-type potentials /16,10/. The collision in- 
tegrals ' ' have been interpolated from the tables of / 37 / (reproduced 
m /10/), yielding the values of r eff and cr eff given in fig.l and Table 1, 
the latter being ~2.4 times as large as the <r eff 's for noble gases. 

Table 1 interaction^ ^ He-H Ne-H Ar-H C-H N-H O-H S-H 
°"ef f (10 cm^) 0.89 1.05 1.45 2.84 2.17 2.35 3.27 

— 7116 neutral ele ment diffusion velocity and time scale. The diffusion ve- 
locity due to the gravitational force for neutral elements is / 6 , 10 / 
v D = D * (A- 1 ) * V(ln p)»r/ n 

v D = 3.70‘1018 .^/t .(A-Dv/TTITa / <r 15 . l/n H • V (In p) • v n (km/day) 

where we have inserted D and n H is in cm -3 and V(ln p) in cm“l. This velocity 
has been evaluated in the framework of model VAL C /43/. r » is close to 1 for 
most high-FIP elements up to x«1900 km (see fig. in / 13/) V except for C (and 
probably S) for which ionization reduces the diffusion velocity by factors up 
to~2.5 within the chromospheric plateau. Fig. 2 shows the calculated v D 's vs. 
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altitude x, for the two extreme elements Ar and C. All other high-FIP elem 
ents have intermediate diffusion velocities. 

From these velocities, we derive for each element A drffusion time 
scales t D A (x) = N H /(S @ v 0fA (x) • n H (x) ) /2/. This time is the time needed for a 
mass of % H- atoms above an altitude x, assumed completely mixed, to be de- 
pleted in atoms A by a factor of 1/e by downward diffusion at its base x 
(S a = solar surface area). These diffusion times t 0 are plotted vs. altitude 
x in fig. 3, for Ar and C. 



5. Discussion - A tentative scheme. The surface of the quiet Sun at altitu- 
des of -700 to 2000 km is believed to be largely covered, above the Super- 
granulation (SG) cells , by horizontal magnetic fields of intensity -5 to 25 G, 
whose energy density exceeds that of the gas above — 1000 km (fig. 4) /3, 11,47, 
14,20,15,22,23/. Above -1000 km, the gas may therefore be stabilized by the 
field, and we may have in the chromospheric plateau a "diffusion region" with 
an horizontal field, which remains stable, isolated from lower lying material, 
over some unknown time scale t . Gravitational settlement of heavy neutrals 
takes place there, but the horizontal field acts as a filter against ions. 
Starting m the upper Chromosphere above the network (SG boundaries) , spicu- 
les are probably the major upward carrier of material from Chromosphere into 
Corona. Almost all the material they lift later falls back, either m the 
spicule itself, or as a ram along the more vertical configuration of the field 
at higher altitudes (fig. 4). Only 0.7 % of it leaves the system m the form 
of SW /above ref./. Above the quiet "diffusion region" at altitude x, we there- 
fore have an apparently fully "mixed region" comprizing the entire TR and Co- 
rona (fig. 4) . The material leaving the "mixed zone" as SW (flux 
treated as constantly replenished by fresh matter from the undelying diffu- 
sion region" at an average effective velocity U(x)= in a 

state fashion lasting over the stability time T of the ^'diffusion region". 

By and large, the abundance of an element A m the "diffusion zone" 

(and hence in the "mixed zone") can be kept low if its downward diffusion 
velocity v D/A (x) relative to the bulk H gas is larger than the effective up- 
ward velocity U (x) of this bulk gas : v D#A (x)/U(x) = n H (x)-v D , A ^ X ^/^SW,H = 

4 > D,A< X >/ 4>SW,A £ 1 (since n A /n H = cst throughout the "mixed zone"). Fig. 2 
shows that this ratio is indeed >1 ( — 30 for Ar, and— 3 for C, our worst 
case) , and does not depend much on the altitude x of the interface between 
the quiet and the mixed-up zones. This strong constraint is thu$ satisfied. 
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m Other constraints lead to limits to the altitude x of the "diffusion 
region. An upper limit x < 2000 kra is imposed by the fact that the high-FIP 
elements must still be mostly neutral /43 , 13/ . A lower limit x £ 1000 km comes 
from the requirement that ions be prevented to diffuse by the horizontal 
magnetic fields ^ addition the stability time of the required 

structure (i.e. the time during which the diffusion region is effectively 
isolated from the denser lower lying material, so that the zones above and 
below it remain unmixed) must be longer than the diffusion time t D (fig. 3) . 
The higher the altitude, the more easily this constraint can be fulfilled. 

For our worst case, C , the required stability time decreases from - 60 days 
at 1000 km to 7 days at 1500 km and a minimum of - 2 days at 2000 km. These 
times are longer than the stability time of the low lying SG and network, 
which is of the order of 1 day /3/. 

6. Conclusion. The heavy element composition anomalies m the TR, Corona, 

SW and SEP can be explained in terms of gravitational settling of the heavy 
neutrals m the Chromosphere, provided that a large fraction of the solar 
surface is covered by quiet, presumably field stabilized, regions at high 
altitude in the chromospheric plateau ( ^ 1500 to 2000 km) • These regions 
must effectively shield the coronal reservoir from the denser, lower chromo- 
spheric gas. They must be stable over at least a few days, and must therefore 
not be strongly perturbed by the daily changes of the SG pattern at lower 
altitude. Unless a different scheme is found to explain it, the odd coronal 
composition therefore permits to set stringent limits to the degree of chaos 
of the Chromosphere, to the turbulent exchange of matter between the lower 
and the upper chromospheric layers . 


Pierre Mem, Georges Michaud and Bernard 
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